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We report electronic transport experiments on a graphene single electron transistor. The device 
consists of a graphene island connected to source and drain electrodes via two narrow graphene 
constrictions. It is electrostatically tunable by three lateral graphene gates and an additional back 
gate. The tunneling coupling is a strongly nonmonotonic function of gate voltage indicating the 
presence of localized states in the barriers. We investigate energy scales for the tunneling gap, the 
resonances in the constrictions and for the Coulomb blockade resonances. From Coulomb diamond 
measurements in different device configurations (i.e. barrier configurations) we extract a charging 
energy of ~ 3.4 meV and estimate a characteristic energy scale for the constriction resonances 
of « 10 meV. 

PACS numbers: 71.10.Pm, 73.21.-b, 81.05.Uw, 81.07.Ta 



The recent discovery of graphene p], 0], filling the 
gap between quasi 1-dimensional (1-D) nanotubes and 
3-D graphite makes truly 2-D crystals accessible and 
links solid state devices to molecular electronics 
Graphene, which exhibits unique electronic properties 
including massless carriers near the Fermi level and po- 
tentially weak spin orbit and hyperfine couplings 0, Q 
has been proposed to be a promisingmaterial for spin 
qubits @ , high mobility electronics 0JH and it may have 
the potential to contribute to the downscaling of state-of- 
the-art silicon technology [9[. The absence of an energy 
gap in 2- D grap hene and phenomena related to Klein 
tunneling [10|, UjJ make it hard to confine carriers electro- 
statically and to control transport on the level of single 
particles. However, by focusing on graphene nanorib- 
bons, which are known to exhibit an effective transport 
gap @, H, EH, [l3[ this limitation can be overcome. It has 
been shown recently that such a transport gap allows to 
fabricate well tunable graphene nanodevices [3, Qj| G3 • 
Here we investigate a fully tunable single electron transis- 
tor (SET) that consists of a width modulated graphene 
structure exhibiting spatially separated transport gaps. 
SETs consist of a conducting island connected by tunnel- 
ing barriers to two conducting leads. Electronic transport 
through the device can be blocked by Coulomb interac- 
tion for temperatures and bias voltages lower than the 
characteristic energy required to add an electron to the 
island [l7[. 

The sample is fabricated based on single-layer graphene 
flakes obtained from mechanical exfoliation of bulk 
graphite. These flakes are deposited on a highly doped 
silicon substrate with a 295 nm silicon oxide layer 
Electron beam (e-beam) lithography is used for pattern- 
ing the isolated graphene flake by subsequent Ar/02 re- 
active ion etching. Finally, an additional e-beam and 
lift-off step is performed to pattern Ti/Au (2 nm/50 nm) 
electrodes. For the detailed fabrication process and the 



sin gle- layer graphene verification we refer to Refs. [14i 
[l8l Il9j|. Fig. la shows a scanning force micrograph of 
the investigated device. Both the metal electrodes and 
the graphene structure are highlighted. In Fig. lb, a 
schematic illustration of the fabricated graphene SET de- 
vice is shown. Source (S) and drain (D) contacts connect 
via 50 nm wide constrictions to the graphene island. The 
two constrictions are separated by w 750 nm and the is- 
land has an area A « 0.06 jim 2 (see Figs. la,b). In order 
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FIG. 1: (color online) (a) Scanning force microscope im- 
age of the investigated graphene single electron transistor 
(SET) device, where the graphene structure and the metal 
electrodes are highlighted. The minimum feature size is ap- 
prox. 50 nm. (b) Schematic illustration of the tunable SET 
device with electrode assignment, (c) Low bias back gate 
trace for Vbi — Vb2 — V pg = V. The resolved transport gap 
separates between hole and electron transport, (d) Effective 
energy band structure of the device as depicted in Fig. lb. 
The tunnel barriers exhibit an effective energy gap of approx. 
6.5 meV. For more information of this model see text. 

to tune the two tunneling barriers and the island elec- 
trostatically and independently, three lateral graphene 
gates [20] have been fabricated closer than 100 nm to the 
active graphene structure (see Fig. la). These are the 
two barrier gates Bl and B2, and the plunger gate PG 
(Fig. lb). The additional highly doped silicon substrate 
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FIG. 2: (color online) Transport as function of the barrier gate potentials Vbi, Vb2 and the back gate at small bias voltages, 
(a) Source-drain current plotted as function of Vbi and Vbi for constant back gate (Vb g = —15 V; see arrow in Fig. lc). Here, both 
individual gaps can clearly be seen. The labels a-c are related to the corresponding close-ups shown in Fig. 3. (b) Symmetric 
barrier gate voltages Vbi = Vbi — 5 V [21[ as function of a varying back gate voltage at V& = 300 fiV. The white areas correspond 
to suppressed current, (c) shows the same for antisymmetric barrier gate voltages Vbi — —Vbi — 5 V, where both transport gaps 
are clearly visible. Please note also the gap homogeneity as function of the back gate, (d) Schematic illustration of the barrier 
configurations explaining the different transport regimes shown in Fig. a. (e) Schematic illustrations for symmetric tuning of 
the tunnel barriers corresponding to Fig. b. (f) Source-Drain current as function of bias and back gate voltage (all other gates 
have been grounded). The measured effective energy gap agrees reasonable well with the model calculation (see arrow). For 
more details see text. 



is used as a back gate (BG) to adjust the overall Fermi 
level (E F ). 

All measurements have been performed in a vari- 
able temperature 4 He cryostat at a base temperature of 
T^ 1.7 K and the sample was heated to 135°C in vac- 
uum for 12 h before cooling down. We have measured the 
two-terminal conductance through the graphene SET de- 
vice by applying a symmetric DC bias voltage V b while 
measuring the current through the SET device with a 
resolution better than 10 fA. For differential conductance 
measurements a small AC bias, V b ^ ac = 50 /iV has been 
superimposed on V b and the differential conductance has 
been measured with lock-in techniques. 

At small bias (V b = 250 /iV < 4k bT) strong current 
suppression is observed at —25 V < Vb g < — 15 V, as 
shown in Fig. lc. This suppression is in agreement with 
earlier studies of graphene nanoconstrictions 0, Q. It 
can be interpreted as a transport gap forming around 
the back gate voltage where the system is charge neutral. 
Hole transport occurs at Vb g < —25 V, electron transport 
at V bg > -15 V. 

Measurements for varying back gate voltage (Fermi 
level) and bias voltage allow to estimate the size of the 
transport gap as shown in Fig. 2f. A value of the order of 
10 meV is found. However, the strong modulation of the 



current shows, that localized states lead to strong trans- 
mission resonances. Therefore we refer in the following 
to an "effective" energy gap or a transport gap. 

The geometric design of our structure (see Fig. la) 
gives local electrostatic access to the constriction regions. 
Fig. 2a shows a measurement of the current where the 
voltages Vbi and V b 2 on the two barrier gates Bl and 
B2 have been independently tuned while the back gate 
voltage was kept fixed at V bg = —15 V. A vertical and 
a horizontal stripe of suppressed current is observed. 
This observation indicates that transport through each 
of the two constrictions is characterized by a transport 
gap which can be individually tuned with the respective 
barrier gate. For example, keeping Vbi = — 20 V constant 
and sweeping Vbi from -20 V to + 5 V keeps constriction 1 
conducting well while constriction 2 is tuned from large 
conductance to very low conductance (into the transport 
gap). The capacitive cross talk from Bl to constriction 2, 
and from B2 to constriction 1 is found to be smaller than 
2%. 

These measurements suggest that the energy diagram 
shown in Fig. Id is a useful description of the data. In 
this figure, high (electron) and low (hole) energy states 
are separated by two solid lines. Outside the constriction 
regions these lines are degenerate and represent the en- 
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ergy of the charge neutrality point in graphene. In the 
constriction regions the two lines are energetically sepa- 
rated indicating the observed effective energy (transport) 
gap Eg by hatched areas. As a result of the lack of an en- 
ergy gap of the two-dimensional graphene material, the 
exact shape of the effective E g (x) (x is the transport di- 
rection) is given only by lateral confinement, i.e., by the 
variation of the width w(x) along the device. We as- 
sume that electron-hole symmetry holds in the confined 
geometry and therefore plot an effective conduction band 
edge at -\-E g (x)/2, and an effective valence band edge at 
-E g {x)/2. 

It is known from earlier experiments 0, @] that 
graphene nanoribbons (or constrictions) exhibit an ef- 
fective energy gap. For ribbons of width w < 20 nm 
the size of this gap scales according to E g = hvp/w, 
where vf = 10 6 m/s is the Fermi velocity. The en- 
ergy gap for nanoribbons wider than 20 nm can be rea- 
sonably well described by E g (w) — a/w ex.p(—bw) [l2| . 
where a = 1 eVxnm and b = 0.023 nm -1 are con- 
stants extracted from fits of the experimental data in 
Ref. [8]. Within this model the width w(x) of our 
graphene structure translates to an effective transport 
band structure exhibiting two tunnel junctions with bar- 
rier height Egfi = 6.5 meV and an almost gap free island 
(E 9i i = 85 /ieV) as shown in Fig. Id. According to 
the model the SET is expected to be operational in the 
regime of \Ep\ < E g ^/2. The measured transport gap 
agrees reasonably well with the modeled barrier height, 
as indicated by the arrow in Fig. 2f. 

The local electrostatic influence of the gate electrodes 
can be incorporated into this heuristic description as a 
local shift of the energy of the charge neutrality point 
described by smooth characteristic potentials <pi(x) (i = 
bl 1 b2 1 pg 1 bg) which may be derived from purely electro- 
static considerations. While 4>b g (x) is independent of x, 
<f>bi(x) and (x) are peaked at the respective constric- 
tions, and (j) pg (x) is peaked within the island. For creat- 
ing the schematic figures in this paper [Figs. 1(d), 2(d), 
(e)] we have used a convenient peaked <j>i{x) function (the 
shape of which is irrelevant for this simple discussion) 
with peak heights compatible with lever arms extracted 
from the experiment (see below). 

Having established a heuristic energy diagram describ- 
ing our sample we now return to the discussion of the 
measurement in Fig. 2a which is facilitated by the di- 
agrams in Fig. 2d. In this measurement Vb g = —15 V. 
From Fig. lc we deduce that the Fermi energy in the con- 
tacts of the structure lies within the conduction band, 
as indicated by the horizontal dashed lines in the four 
drawings in Fig. 2d. The four drawings represent energy 
diagrams corresponding to the four corners of Fig. 2a as 
indicated by the white numbers. In corner 2 transport 
takes place in the conduction band throughout the whole 
structure. In corner 1 (4) transport occurs in the conduc- 
tion band in the right (left) part of the structure. The left 
(right) constriction is traversed via states in the valence 
band. The situation is even more complex in corner 3, 



where the Fermi energy cuts both barrier regions in the 
valence band. Although these situations imply two or 
even four p-n-like transitions along the structure, no dis- 
tinctive features are observed in our measurements. This 
may be a manifestation of the suppression of backscat- 
tering due to Klein tunneling. 

Figs. 2b and 2c demonstrate the consistency of our 
heuristic model with the experimental observations. 
Fig. 2b shows the current measured as a function of Vb g 
and Vbi, with Vb2 being simultaneously swept such that 
Vb2 = Vbi + 5V (see dashed line in Fig. 2a). In this 
way the barrier regions are simultaneously shifted up or 
down (see Fig. 2e). Fig. 2b shows that the transport 
gap measured as a function of the back gate is shifted 
correspondingly, with AV^/AV&i^ ~ 0.9. 

Fig. 2c shows the current measured as a function of Vb g 
and V51, with V52 being simultaneously swept such that 
V b i + V b2 = 5 V (see dotted line in Fig. 2a). For V bl = 
±15 V (vertical dashed lines in Fig. 2c) the position of the 
gaps in energy correspond to diagrams 1 and 4 in Fig. 2d. 
In these two cases, sweeping the back gate allows to probe 
the two spatially separated transport gaps individually. 

If we focus on a smaller voltage scale much more 
finestructure in the Vbi — Vbi parameter plane appears, 
as shown in Fig. 3. Subplots 3a-c are different close ups 
of Fig. 2a (see black labeled boxes therein). Although 
Figs. 3a-c show the current in three different regimes the 
transport characteristics do not differ significantly. Here, 
we distinguish between the PP (Fig. 3a), NN (Fig. 3b) 
and the NP (Fig. 3c) regime, depending on either hav- 
ing the tunnel barriers (according to Bl and B2) shifted 
down (N) or up (P). We observe in all regimes (Fig. 3) 
sequences of horizontal and vertical stripes of suppressed 
current and current resonances. Their direction in the 
Vb\-Vb2 plane indicates that their physical origin has to 
be found within constriction 1 (vertical stripes) or con- 
striction 2 (horizontal stripes). A blow-up of a small 
region in Fig. 3a is shown in Fig. 3d. The current ex- 
hibits even finer resonances which are almost equally well 
tuned by both constriction gates. We therefore attribute 
these resonances to states localized on the island between 
the barriers. It will be shown below that these reso- 
nances occur in the Coulomb blockade regime of the is- 
land. We attribute the deviations from perfectly straight 
diagonal lines to the presence of rough edges and inho- 
mogeneities within the graphene island which has dimen- 
sions (slightly) larger than the elastic mean free path. 

This characteristic pattern (Fig. 3d) can be found 
within a large Vb\-Vb2 parameter range within the regime 
where the two barrier gaps cross each other (i.e. the inner 
bright part of Fig. 2a). 

So far we mainly focused on the barriers and in the 
following we concentrate on the charging of the island 
itself. We fix the barrier gate potentials (Vbi and V b 2) 
either in the NN regime or in the NP regime in order to 
study Coulomb blockade. Fig. 4a shows sharp conduc- 
tance resonances with a characteristic period of about 
20 mV (Hi =5.570 V, and H2 =-2.033 V are fixed). Their 
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FIG. 3: (color online) Source-drain current through the 
graphene SET as function of the barrier gates Vbi and Vhi 



for constant bias Vb = 300 jiY and back gate Vb g 



-15 V. 



(a-c) are close-ups of Fig. 2a (as indicated therein by labeled 
boxes), showing transport in the PP (a), NN (b) and NP (c) 
regime. Ontop the horizontal and vertical transmission mod- 
ulations we observe (diagonal) Coulomb blockade resonances. 
This is best seen in Fig. d, which is a close-up of Fig. a. In 
Figs, a and d the current has been multiplied by a factor 2 
and 10, respectively, to meet the color scale shown above b. 



amplitude is modulated on a much larger voltage scale of 
about 200 mV by the transparency modulations of the 
constrictions (cf. Fig. 3d). These resonances in the nar- 
row graphene constrictions can significantly elevate the 
background of the Coulomb peaks (see e.g. black ar- 
row). The inset of Fig. 4a confirms that transport can 
also be completely pinched off between Coulomb block- 
ade peaks. Corresponding Coulomb diamond measure- 
ments [I?}, i-e., measurements of the differential conduc- 
tance (Gdiff = dl/dVb) as function of bias voltage Vb and 
plunger gate voltage V pg are shown in Fig. 4b. Within 
the swept plunger gate voltage range no charge rearrange- 
ments have been observed and the peak positions were 
stable over more than 10 consecutive plunger gate sweeps. 

In Fig. 4c we show conductance resonances, which 
have been measured within the NN regime (for fixed 
V bl = 8.79 V and V b2 = 8.85 V, see Fig. 3b). The V pg 
range shown here is wider than in Fig. 4a. Again we ob- 
serve (i) strong transport modulations on a V pg scale of 
about 100 mV, which originate from resonances within 
the barriers and (ii) Coulomb peaks on a V pg scale of 
about 20 mV which are blown up in Fig. 4d. The corre- 
sponding Coulomb diamond measurements (Fig. 4e) are 
similar to those measured in the NP regime (Fig. 4a). 
The Coulomb peaks (Fig. 4d and inset in Fig. 4a) and 
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FIG. 4: (a) Source drain current as function of the plunger 
gate voltage V pg at fixed back gate and barrier gates in the NP 
regime (V bg = -15 V, V bgl = 5.67 V, and V bg2 = -2.033 V). 
The inset (close up) clearly shows Coulomb peaks (b) Cor- 
responding Coulomb diamonds in differential conductance 
Gdiff , represented in a logarithmic color scale plot (dark re- 
gions represent low conductance). A DC bias Vbias with a 
small AC modulation (50 jiV) is applied symmetrically across 
the dot and the current through the dot is measured, (c) 
Coulomb resonances ontop and nearby strong transport mod- 
ulations in the NN regime (V bg = -15 V, V bg i = 8.79 V, and 
Vb 9 2 = 8.85 V). (d) shows a close up highlighting Coulomb 
peaks and (e) shows the corresponding Coulomb diamond 
measurements. The color scale is adapted from Fig. b. 



the Coulomb diamonds are not very sensitive to the tun- 
nel barrier regime, although in one case a p-n-like junc- 
tion should be present, whereas in the other case a more 
uniform island is expected. 

From the extent of all the diamonds in bias direction 
we estimate the average charging energy of the graphene 
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TABLE I: Capacitances and lever arms of the different gate 
electrodes, including source and drain contacts, with respect 
to the graphene island. Most values are independent from the 
measurement regime, NN or NP. If there is a difference the 
NP value is given and the NN value is put in brackets. 
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FIG. 5: Nearest neighbor Coulomb peak spacing statistics in 
both the NN (left) and NP (right) regime. Several plunger 
gate sweeps (at different barrier gate potentials) have been 
evaluated and in total 431 Coulomb peaks in the NN (and 
460 peaks in the NP) regime have been considered for the 
statistics. 



single electron transistor operated in both regimes to 
be Ec ~ 3.4 meV. This charging energy corresponds 
to a sum-capacitance of the graphene island = 



1 /E C 



47.3 aF, whereas the extracted back gate 



capacitance C b g ~ 18 aF is higher than the purely geo- 
metrical parallel plate capacitance of the graphene island 
C = eoeA/d w 7.4 aF. This is related to the fact that the 
diameter of the graphene island (y/A) is approximately 
the same as the gate oxide thickness d [HI, [22[ . 

The lever arms, and the electrostatic couplings of the 
electrodes to the graphene island do not change signif- 
icantly between the NN, PP (not shown) and the NP 
regime. Thus, the lever arm of the plunger gate is 



C P g/Cj: « 0.15 (C p 



6.9 aF), whereas the 



a P9 

electrostatic coupling to the other gates were determined 
to be C bl « 5.5 - 6.0 aF and C b2 ~ 5.0 aF. All 
lever arms and capacitances are summarized in Tab. 1. 
It shows that the island geometry and dot location with 
respect to the lateral gates stays almost constant. How- 
ever, the capacitive coupling to the source and drain con- 
tacts (i.e. Cs and Cd) changes significantly as function 
of the tunnel barrier configuration. This can be nicely 
seen when comparing the symmetry of the diamonds in 
the NN and NP regime as shown in Figs. 4e and 4b. 
While the size and fluctuations of the diamonds remain 
(almost) constant the lever arms of the source and drain 
contacts change strength. In one case (NP regime) we 
extract Cs ~ 1.8 aF and Cd ~ 9.6 aF, whereas in the 
other (NN regime) C s ~ 10.1 aF and C D « 1.8 aF, 
which can be seen from the different slopes of the di- 



amond edges. However, the individual tunnel barriers 
strongly depend on the local barrier configuration and 
change also within the NN or the NP region. 

We now estimate the energy scale of the resonances 
in the constrictions. The spacing of the constriction res- 
onances in plunger gate is about 200 mV, whereas the 
spacing of Coulomb peaks is 20 mV. By assuming that 
the capacitance between the plunger gate and the local- 
ized states in the constrictions leading to the resonances 
is about three times smaller than C pg (estimated from 
the geometry of the device) the energy scale of the reso- 
nances in the constriction is about 10 mV, in agreement 
with the measured gap in Fig. 2f. 

Alternatively, this characteristic energy scale can also 
be estimated by considering that the back gate voltage 
sweep from -25 V to -15 V (around the charge neutrality 
point at Vh g = -20 V, Fig. lc) translates to a Fermi en- 
ergy sweep over an energy interval of approx. 120 meV. 
Near the Dirac point the spacing of the constriction res- 
onances in back gate voltage is found to be of the order 
of 200 mV leading again to a characteristic energy scale 
of 10 meV. 

Finally, we also performed Coulomb peak spacing 
(AV) statistics in both, the NN and NP regime with in 
total more than 900 Coulomb peaks, as shown in Fig. 5. 
The mean nearest neighbor spacing of the Coulomb peaks 
in both the NN and NP regime do not differ significantly 
(AFatat=17.4 mV and AFatp=17.9 mV). The broad- 
ening of the peak spacing distribution is in both cases 
significant ranging from <jnn ~ 3.3mV (0.6 meV) to 
(Jjsfp ~ 2.5mV (0.5 meV), which is in agreement with 
Ref. UK . The broadening of the observed unimodal peak 
spacing distribution is significantly larger than the dif- 
ference between the average spacings in the NN and NP 
regimes (0.075 meV). The inhomogeneity of the island, 
as indicated by the different slopes in Fig. 3d may sig- 
nificantly contribute to the observed broadening, which 
might be also partly influenced by the underlying modu- 
lation of the transmission through the narrow graphene 
constrictions. The broadening of the distributions is 
significantly larger than expected for a purely metallic 
SET [23[ . On the other hand the width of the distribution 
is of the order of the estimated single-particle level spac- 
ing similar to previous observations in high-quality 
GaAs quantum dots [25|, [26[ . This may indicate the im- 
portance of quantization effects. 

In conclusion, we have fabricated and characterized a 
fully tunable graphene single electron transistor based 
on an etched width-modulated graphene nanostructure 
with lateral graphene gates. Its functionality was demon- 
strated by observing electrostatic control over the tunnel- 
ing barriers. From Coulomb diamond measurements it 
was estimated that the charging energy of the graphene 
island is « 3.4 meV, compatible with its lithographic 
dimensions. These results give detailed insights into tun- 
able graphene quantum dot devices and open the way to 
study graphene quantum dots with smaller dimensions 
and at lower temperatures. 



6 



Acknowledgment — The authors wish to thank sions. Support by the ETH FIRST Lab and financial 
R. Leturcq, P. Studerus, C. Barengo, P. Strasser, support by the Swiss National Science Foundation and 
A. Castro-Neto and K. S. Novoselov for helpful discus- NCCR nanoscience are gratefully acknowledged. 



[1] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 
M. I. Katsnelson, S. V. Dubonos, I. V. Grigorieva, 

A. A. Firsov, Science, 306, 666, (2004). 

[2] For review see: A. K. Geim and K. S. Novoselov, Nat. 
Mater. 6, 183 (2007), A. H. Castro Neto, F. Guinea, 
N. M. Peres, A. K. Geim, cond-mat 0709.1163vl (2007). 

[3] C. Joachim, J. K. Gimzewski and A. Aviram, Na- 
ture, 408, 541, (2000). 

[4] H. Min, J. E. Hill, N. A. Sinitsyn, B. R. Sahu, L. Klein- 
man, and A. H. MacDonald, Phys. Rev. B, 74, 165310, 
(2006). 

[5] N. Tombros, C. Jozsa, M. Popinciuc, H. T. Jonkman and 

B. J. van Wees, Nature, 448, 571-574, (2007). 

[6] B. Trauzettel, D.V. Bulaev, D. Loss, and G. Burkard, 

Nature Physics, 3, 192, (2007). 
[7] Z. Chen, Y. Lin, M. Rooks and P. Avouris, Physica E, 40, 

228, (2007). 

[8] M. Y. Han, B. Ozyilmaz, Y. Zhang, and P. Kim, Phys. 

Rev. Lett., 98, 206805 (2007). 
[9] M. Ieong, B. Doris, K. Kedzierski, K. Rim and M. Yang, 

Science, 306, 2057 (2004). 
[10] N. Dombay, and A. Calogeracos, Phys. Rep., 315, 4158 

(1999). 

[11] M. I. Katsnelson, K. S. Novoselov, amd A. K. Geim, Na- 
ture Phys. 2, 620625 (2006). 

[12] F. Sols, F. Guinea and A. H. Castro Neto, Phys. Rev. 
Lett., 99, 166803 (2007). 

[13] X. Li, X. Wang, L. Zhang, S. Lee, H. Dai, Science, 319, 
1229 (2008). 

[14] C. Stampfer, J. Giittinger, F. Molitor, D. Graf, T. Ihn, 
and K. Ensslin, Appl. Phys. Lett., 92, 012102 (2008). 

[15] L. A. Ponomarenko, F. Schedin, M. I. Katsnelson, 
R. Yang, E. H. Hill, K. S. Novoselov, A. K. Geim, Sci- 
ence, 320, 356 (2008). 



[16] F. Miao, S. Wijeratne, Y. Zhang, U. C. Coskun, W. Bao, 

C. N. Lau, Science, 317, 1530 (2007). 
[17] L. P. Kouwenhoven et al. 1997, "Electron transport 

in quantum dots", Mesoscopic Electron Transport (ed) 

L. L. Sohn, L. P. Kouwenhoven and G. Schon (NATO 

Series, Kluwer, Dordrecht). 
[18] A. C. Ferrari, J. C. Meyer, V. Scardaci, C. Casiraghi, M. 

Lazzeri, F. Mauri, S. Piscanec, D. Jiang, K. S. Novoselov, 

S. Roth, and A. K. Geim, Phys. Rev. Lett. 97, 187401 

(2006). 

[19] D. Graf, F. Molitor, K. Ensslin, C. Stampfer, A. Jungen, 
C. Hierold, and L. Wirtz, Nano Lett., 7, 238 (2007). 

[20] F. Molitor, J. Giittinger, C. Stampfer, D. Graf, T. Ihn, 
and K. Ensslin, Phys. Rev. B 76, 245426 (2007). 

[21] Here an offset of 5 V has been used in order to com- 
pensate the slightly different doping of the two narrow 
graphene constrictions. This offset can also be nicely seen 
at Vb2 in Fig. 2c. 

[22] See Fig. 12.11 in T. Ihn, Springer Tracts in Modern 
Physics 192, Springer 2004, p. 102. 

[23] M. Furlan, T. Heinzel, B. Jeanneret, S. V. Lotkhov and 
K. Ensslin, Europhys. Lett., 49, 369 (2000). 

[24] The single level spacing is estimated by A qm = 
Hvf /2y AN/tt, where A is the area of the graphene is- 
land and N is the number of carriers on the island. For 
N=10 to 100 electrons on the island we estimate a char- 
acteristic energy spacing of 0.75 meV to 0.24 meV. 

[25] S. R. Patel, S. M. Cronenwett, D. R. Stewart, A. G. 
Huibers, C. M. Marcus, C. I. Duruoz and J. S. Harris, Jr., 
K. Campman and A. C. Gossard, Phys. Rev. Lett. 80, 
4522 - 4525 (1998) 

[26] S. Liischer, T. Heinzel, K. Ensslin, W. Wegscheider, and 
M. Bichler, Phys. Rev. Lett. 86, 2118 (2001) 



